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A MEAN-VALUE THEOREM OF THE 
RIEMANN ZETA FUNCTION 


By F. T. WANG (Hangchow) 
[Received 30 July 1945] 
1. Ir is well known that, for any fixed o > }, 


T 
z{ IS(o-Fit)|? dt ~ (20) (T+ 0). 


[ shall establish in this note the following analogous theorem: 
T 
1 , ' 
z/ |f(o+-2t) |? log |f(o+-it)| dt ~ f(2c)log £(20). 


0 
Proof. By Jensen’s inequality,* when we take 


P= (Cl, f= (cl, P(x) = rlogz, 
T 
T \C|? dt 
we get | \£|*logiz| de > if If? dtlog ® —__, 
0 \f|2-€ dé 
si 
= k—1\2 ) 
Let us put C,(¢) = sit p ori. 2 
et us pu Ko) IT{>( : ya j (2) 
By a theorem due to Titchmarsh and Ingham, fer an arbitrary given 
positive fixed number « we can determine a number 7. such that 


T 
1 2 
T | |g|?-€ dt —C,_,.(20)| < €* 
0 
when 7’ > T.. — 


fige- dt = TC, (20) +O(Te) 
0 
when 7’ > 7. From (1) and (3) we have 
T 
J iglttogig| ae > +{TE(20 )+O(eT)}log Rodan x *) 
0 C,- 3e(20 a)+ Ole ¥ 
¢(2c) 
> TC(2 oi O 4 
> TE(20) log a= 95) + Oe) (4) 
when 7 > 7. 
* Hardy, Littlewood, and Pélya (1), 151. + Titchmarsh (3), Ingham (2). 
3695.18 B 
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e\? e \? 

‘i (\-: \0-9) 1-3) 
4 2v 
1—eL(v)+ O(e?2” log v), 


where 


Therefore 


log 


@ 
since > > p*-* log v is convergent. 
pv=1 


By (2) and (5) we have 


log aC eo? log > 


1-te(2 


= e€log {(20)+- O(e?). 


Hence, when 7' > T,, we have 
- 
[ \Z\Plog|g| dt > TL(2e)log £(20)+O(Te). 
f 

Taking the successive limits 7’ > oo and « > 0, we have 


. 
lim 5 | ||? log || dt > {(2e)log f(2c). (7) 


To 


I now use a general inequality concerning the arithmetical and 
geometrical mean,* in which I take 


p=(t?, f= ltl r=e. 


few a fe dt ; 
__—______ < log{° — ; 
sd ight 
fi |? de f eat 


* Hardy, Littlewood, and Pélya (1), Th. 184. 
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T 
C1 44(2¢) 
Hence J \¢|Plog|<| dt < T(20)— liog oe 5 + 0(Te) 
when 7 > T.. 


As in (6), we similarly get 
T 
J [Sllog|g| dt < TE(2e)log £(20)-++ O( Te) 
a 
when 7 > 7. Hence 


T 
Tim 7, | IL|Plog|£| dt < £(20)log £(20). (8) 
0 


The theorem follows from (7) and (8). 
2. Hardy and Littlewood* proved that, when o = }, 


T 
f SG+it)? dt ~ Tlog 7. 
0 


I cannot extend my method to establish the analogue. The follow- 
ing inequality, however, is true for sufficiently large values of 7’: 
T 
A,T log T log, T < f |€(4+-%t) |? log |f(4-+-2t)| dt < A, T log T log, T. 
0 
In general, we can prove that, when 0 < k < 2 and oa > 3, 





— k—1 
r 2 Gir )Bo-% 
[ \S(o-+it)|*logiZ(o+it)| dt ~ TO,(20) 5 =, ; 
a Dp >: v+k—1\? —2vo 
v=0 z ‘ 
where E(v) = ee 
r= a 


* Titchmarsh (4). 


REFERENCES 
Hardy, Littlewood, and Pélya, Inequalities (Cambridge, 1924). 
. A. E. Ingham, ‘Mean-value theorems and the Riemann zeta function’: 
Quart. J. of Math. (Oxford), 4 (1933), 278-90. 
3. E. C. Titchmarsh, ‘Mean-value theorems in the theory of the Riemann 
zeta function’: Messenger of Math. 58 (1929), 125-9. 
The Zeta Function of Riemann (Cambridge, 1932). 
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ON THE ZEROS OF THE RIEMANN 
ZETA FUNCTION 


By E. C. TITCHMARSH (Ozford) 
[Received 18’ June 1946] 


1. Let N(7’) denote the number of zeros B+iy of {(s) such that 
0< ¥< < 7, and let N,(7’) denote the number of zeros such that 
B=}3,0<y<T. The Riemann hypothesis is equivalent to the 
statement that N,(7') = N(7) for all values of 7’. 

It is known that N(7’) ~ (27)T log T as T’' +> oo. It was proved 
by Hardy and Littlewood* that 


N,(T) > AT (1.1) 
(A denotes various absolute constants). Recently, A. Selbergy has 
shown that N(7) > AT log T. (1.2) 


This is a remarkable improvement, since it shows that a finite propor- 
tion of the zeros of f(s) lie on the line R(s) = 3. Actually Selberg also 
propes Shes Nj(T+U)—N(T) > AU log T (1.3) 
if U > T*, where a > }. 

The object of the present paper is to give a shorter proof of (1.2). 
This is done by using a Fourier transform method similar to that 
which I used in my Cambridge Tract to prove (1.1). Apart from this, 
the essential ideas are of course taken from Selberg’s paper. 

I communicated my proof privately to Dr. Selberg, and in reply 
he pointed out that a considerable part of it was unnecessary for the 
proof of (1.2), since it followed the lines of his more detailed analysis 
used in proving (1.3). The modification due to Dr. Selberg is incor- 
porated in this paper, the analysis suggested by him beginning at § 9 

2. The method of Selberg consists of modifying the series for C(s) 
by multiplying it by the square of a partial sum of the series for 
{{(s)}-*. We define a, and B,, by 


er: (@=oc+it;>o>1), y=], 


* G.H. Hardy and J. * Littlewood, ‘The Zeros of Riemann’s Zeta Function 
on the Critical Line’, Math. Zeitschrift, 10 (1921), 283-317. See also my 
Cambridge Tract, The Zeta Function of Riemann, Theorem 32. 

+ A. Selberg, ‘On the Zeros of Riemann’s Zeta Function’, Norske Viden- 
skaps-Akad. Oslo, Mat.-Naturv. Klasse, 1942, No. 10, 1-59. 
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and i= o,(1— 7087) (l<v< &). 


Since a IT ('-3) a Md (orb 


we have Oy, Ay = ay », 


if (v,,¥2) = 1. Since the series for (1—z)* is majorized by that for 
(1—z)-*, we see that, if 


then |«,| < a, for all values of v. In particular |«,| < 1, |B,] < 1 for 
all v. 


3. Let d(s) = 2,» , 


where € is to be determined later as a function of 6 (or of 7’). Let 
e+e 
I'(48)a-#8f(s)h(s)6(1—s)z8ds (c > 1). 
c—iw 
There is a pole at s = 1 with residue z¢(1)4(0). Hence 
$+ico 
O() = AHO) +E; [Pee HC(2)4(6)6—e)e de 


$i 


Now, if s = 4+4#, 
T'(}8)n-#¢(s) = —25(t)/(P+}), 


where 3(¢) is real for real ¢. Hence 


i) 


(2) = e690) — 2 i | E00 ipa initet at 


On the other hand 


I\(48)z —}s re 
neusy -s 
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where p and » (and later x and A) always run from 1 to é. Hence 


Je © 1 4(4-+-ét)| 22H at 


+ 
= 4244(1)6(0)— > vexp( “EF. 


Putting z = exp{—i(}n—}5)—y}, it ‘follows that the functions 
1 X(t) 
(27) P+-4 


fly) = ¥249()$0)—24 YY BeBe exp(— 7) 
n=1 mh Vv 


are Fourier transforms. Since F(t) is real, it follows that 


F(t) = $($-+it)| 2e't7 480, 


o +H 


{tf F(u yan) ae = = sf vw sm ay, (8.1) 


—o t 
where H is to be chosen later (it is to be small, and it may be assumed 
that it is less than 1). The right-hand side is less than 

1/H 


2H? J If? dy +8 i If y)l?y- dy. (3.2) 


1/H 
Putting y = logx, G = e¥4, the first integral in (3.2) is equal to 


—i(jn— re 0-> > he BB, rexp(—™ 7m ~ TL" pitt —8)72 ‘ dz, 
Oe 


— the triple sum ne ), this is not greater than 


G 
sf ae ao ~ de { ae < ssa? | lg(x)|* de. 


i 
Similarly the second integral in (3. 2) does not exceed 


gl F Ig(z)l? 


} log’a 


We have to obtain upper bounds for these integrals as § > 0, but 
it is more convenient to consider directly the integral 


J(X,0) = | lg(x)|22-9 da (0<0<4; X>1). 
x 
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4. We have 
J(X, 6) 


S > ee am 
ZEEE oleae 


m=1 n=1 
2,,2 
(4.1) 


Let >, denote the sum of those terms in which mx/A = nyu/v, and >, 
the remainder. Let (xv,Au) = q, so that 

Kv = aq, Au = bq, (a,b) = 
Then, in >, ma = nb, so that n = ra, m = rb (r = 1,2,...). Hence 


x reiki)” 
The last r-sum is of the form 


aa -jexm 


where K(6), and later K,(@), are bounded functions of 6. Hence we 
obtain 


r 0-a[ ¢ 
ama ew dy +012) | | ev'y~ dy +0fc9)-4] 4 
0 0 


2] 


re K(0) [ e-vty-6 dy +0(.vn] + O{X1- log(Xvn)} 


0 





Ny K,(6)n9-+ Xxi-0 
= ~xnit— on +0{=—log(Xvn)}. 


Putting 7 = 27xu®q-*sin 5, it follows that 


= S(0) K,(@) 
21 = a@sinsyiaxot — 9 


1-0 Av: - vl 
fe > 3 > > ip Pan P| (4.2) 


(2a sin §)#9-#S(@) + 
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where S(é) = ‘ 2 > ps (z)" "BePababe, (4.3) 


Let ¢_o(p) = p'~? me = pi “TT (-; =): 


m\p pip 


where y(m) is the function of Mébius. Then 


qi = 2, $-0(e) = oi r o(p >. ¢-0() 


” plas plAw 
Hence 


0) = F 4-0) D> wir >to oe) 
(4.4) 


Let d and d, denote positive integers whose prime factors divide p. 
Let « = dx’, v = d,v’, where (x’, p) = 1, (v’,p) = 1. Then 


2 Eat = oi > Bae, > = Bay. (4.5) 


pl«v pldd; w'<é/d wv’ <b/d; 
Now for («x’,p) = 1 
J log{é/(dx')} a 
Bax’ =e log € _ ~ Jog é° et 


Hence (4.5) is equal to 


a a ] 3 oy] g 4.6 
oats 2¢ z d'- a 1-6 Fa, a v’ a omy 


plddy Kn’<é/d v'<él 





5. We now prove three lemmas. 
LEMMA a. We have 
0? e g 5 
> Ap 8 F5 = of (6 ‘) log? = +11 (1+ | (5.1) 
K’<&/d p\p 
uniformly with respect to 0. 


We may suppose that € > 2d, since otherwise the lemma is trivial. 


Now 
1+io : 
J © ds = § (0<#<}), 
Qart 3? \loga («> 1). 
1—io 
Also 
On 1 \3 ( 1 te l 
ee 1 oo eG = 1 Dae pee, SE - 
tle TL per) I] pits] {(1—8+8) 








\ 
H 
| 
H 
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Hence the left-hand side of (5.1) is equal to 


STW 4 


to pip 





There are singularities at s = 0 and s = 0. If # > {log(é/d)}-", we can 
take the line of integration through s = 6, the integral round a small 
indentation tending to zero. Now 

1 
C(1+1t) 
for all ¢ (large or small). Also 


T(t 5)~ LT +2) 


Pip pip 


< Alt 








}= ALT +3 


pip 








Hence (5.2) is 
(2°17 [43.28 fa (e\? a) 
& Oe Been: at om £ Pe 
ol(a) IT +5) e+e ol(a) II +5) Br” 
p\p a p\p 
and the result stated follows. 
If 6 < {log(é/d)}-, we take the same contour as before modified by 


a detour round the right-hand side of the circle |s| = 2{log(€/d)}-. 
On this circle 


\(€/d)*| < & 
i. 4 1\3 
rT 5-0) : oT T(:+5) | 
as before, and \¢(1—6+8s)| > A log(é/d) 


since &() has a simple pole at w = 1. Hence the integral round the 


© fee sSTT (432 fle) ~ fue TT (+3) 


pip p\p 








The integral along the part of the line o = @ above the circle is 


g\8 lk of dt é 

A TO+A I a)- Glee Tp 
Pip A(log €/d)>* 

The lemma is thus proved in all cases. 


Lzemoa f. Para = Ol. T] (1 +°)|. 


piddy pip 
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Defining «jz as in § 2, 
>, 33 nee oa %a,| Oa Xd, y 1 
plddy <2 dd, pilD D 
(where D is a number of the same class as d or d,) 


= TI0-3)"= (3) 


Pip Pip 
the result stated. 


LEMMA y. S(@) = o( i) 


uniformly with respect to 6. In particular 


1 
810) = O(ez) 
By (4.5), (4.6), and the above lemmas 


SSM oh 5 el me fee TT 
- SeT Pa 
“= PT (+ *) | 

= = Ol i TI (+5) ] 























p<é pip 
= (+5) 
= i+ 
Ol eae SP a II 
£20 1 1 
= Olea f =) 
. log*é the pit <- ni 
since 
1\4 4 1 1 
[7 bey)- oll bs) Arnal a) 
pip pip pip n|p 
Hence 
£20 1 £0 Ss 1 l ! 
S() = Ol = oe ey 4 
( ) log’é &y, othe (mp,)**nt logé = ni+é ake pito 





> at 7 lege) 


pe P 


= of Se 23 ni 
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6. By Lemma y, the first term on the right-hand side of (4.2) is 


1 
Ol soxe log : (6.1) 
The second term is 


CF igs) ~ °(soxtog) 


if XE? = O(8-+). The third term is 
of 7 > > > ; : 7} = ofiales log(XE/8) a1, e%¢| 


1 
sf eax a 
1 
7" Ds _ 
provided that Xé 5 ata (6.2) 
We shall ultimately choose é so that (6.2) is satisfied. It will then 


follow from (4.2) that 
1 
2d certs sai 


7. Now consider >,. If P and Q are positive, and X > 1 


r dx 1 f e-Pu : e-P 
—Pz*+iQz* sit aad be — 2 
| . ~~ 3 jini x) 
x? 


e.g. by applying the second mean-value theorem to the real and 
imaginary parts. Hence 


%- lpr LLIB LD Pe as "El x 


co-op 


The terms with mx/A > nu/v contribute to the (m,n) sum 
of S e—7m*x*)-* sind zz r-") | 
? v2 : 
m=1 n<m«v/Au 


A v 


| ee ali 


m2 nu? . mK “(m mH) mx(mxv—NAu) 
a 


A*v 
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and > . : 1+ — at BR aun = = 140°"), 


— mK — mi 


Hence the (m,7) sum is 


olke S (+ mee arm? K2A- = 


K m mA 
\ m=1 id 


A? meth x 1 2 
1 = lnot 
mi © {( "; ‘log ia ind +h, 8 and | 


3 a ] 
log =-+ — log? = 
(08 +5, ” 3) 
provided that log € = O(log 1/5). The terms with mx«/A < np/v may 
be treated similarly. It follows that 


= wA o* | 
3.= ogo D DD (owg+Liry) 
== o(x-% log é log 5 + X-9€? log’é log? 3 


= O(X-*étlog* 5), (7.1) 


8. Lemma 8. If & = 8-, where c is a constant less than }, and (6.2) 
is satisfied with X = G, then 


t+H 


r H 
Jl, J Fe u) du "aa Ol soz (8.1) 


From (6.3) and (7.1) it follows that 
asec dam Sy 


uniformly with respect to 6. Hence 


fice \g(x == f+ fae = = [-—2 J]; 40 fa LJ dx 


a a _ af logG 
- ew { watt) -_ (5 e*). 
1 


taking e.g. 6 = 4. Also 
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i 6J(G, 0) dd = J Ig(a)|2 da j 6-9 do 





fi 1 tv 
ge Geers ~ Qatloga 2 log®x 


ste)" 
d: 
ioe -[ 35 


if G is sufficiently lines. a 


3 
dé 1 
Ee oa 
0 


1 
ols log G log } 
Also 4(0) = O(€), d(1) = O(log €). The result therefore follows from 
the formulae of § 3. 


9. So far the integrals considered have involved F(t). We now turn 
to integrals involving |F(é)|. The results about such integrals are 
expressed in the following lemmas. 


LEMMA e. f |F(t)|2 dt = o(5; sora 


By the Fourier-transform formulae, the left-hand side is equal to 


r r e—ita—48) 2 
2 { ifanitay = 2 [ $(1)4(0)—g(a)} de 


2x 
0 


<4] |g(x)|? dx + O(€* log*é) 


in the notation of § 3. 
The integral here is J(1,0). This can be estimated in much the 
same way as J(X,@). We follow the analysis of § 4 until we come to 


the sum 1 7 
>: — = log $4740 “1 
r 


rsylvn 
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This leads to a sum similar to S(0), but containing a factor log(ky/q). 
In dealing with the logg, we have to replace the previous ¢_(p) by 


| —set-ole)| = (loge + > 28%) (p). 


pip 
The factor in brackets is O(log p) = O(log 1/5). Terms of the same 
order arise from log K and logy. The result is therefore similar to 
(8.2), but with 1/(@X®) replaced by log 1/8. 
This proves the lemma. 


t+H 


H? log 8-1 
i in = a. ee 
LEMMA ¢ siti |F(u ( o(* Slogé 
For the left-hand side does not exceed 
S t+H 
it yh w|? du} de H | FCW ‘aw | a= HE f |F(w) )|2 du, 


a —@ 


and the result follows from Lemma e. 


Lemma 7. [fd = 1/7, 


: 
J |F(t)| dt > AT?. 
0 


2+4 2+i7 3+iT +7 
We have | + J + J + | }etergre) ds = 
b+i ott HiT BIT 
Now ¢(s) = O(&!) for o > 4; hence the first term is O(€), and the 
third is O(€7'*). Also 
— a 
c(s)p%(s) = 1+ > =, 
n=2 
where |a,,| < d;(n), the number of ways in which can be expressed 
as the product of three factors. Hence 





2+i/' _ 2+iT 
ds 
J Here) do =i T—1)+ > a Ie 
i i(7—1)+0 3 ten 
=i7T+00). 


- y 
It follows that | L(4+- it) 6?(4-+ it) dt = T+0(7). 
0 
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Hence 
- - 
[ FO) at > A f -1L4+ingr(gt in) de 
0 0 
T 
> AT [ |C(b-+it)g2(4+it)| dt 
+7 
¥ 
> ATA) [ ((4-+it)g(h+it) det 
4T 
> AF. 
i t+H 
Lema 0. [at | |F(u)| du > AHT?. 
0 t 


The left-hand side is equal to 


T+H min(7',u) ¥ u y 
[ |F(u)| du f dt > | |F(u)| du f dt = H f |F(u)| du, 
0 max(0, u—H) H u—H H 


and the result follows from Lemma 7. 


10. We can now proceed to the proof of (1.2). Let Z be the sub-set 


of (0, 7’) where 
t+H 


J |F(u)| du > 
t 





[f this holds, then F(u) must change sign in (¢,¢+-H); hence so must 
E(w), and so ¢(4-++iw) must have a zero in this interval. Now the 
above inequality gives 

t+H t+H 


jr 


t 











| 








The left-hand side is not greater than 


if af (ic w)| du) ail? < < Vm(E By f (rs u)| au) ar\* 


< A{m(E)}\HT* 
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by Lemma ¢, with 6 = 1/7, £ = T°. The second term on the right 
is not greater than 


[fas] [#0 


o,t+H 


du| at}* < < mf J F(u) du 


AT*H* _ ATH 
Stlogté ~— logté * 


vat}? 














Hence Vm(E£) > A, Tt—. 4.7m : 
H* log*é 
where A, and A, denote the particular constants which occur. Taking 
44 4 
Ajlogé AjclogT 
it follows that m(EH) > tAeT7. 


This choice of H implies that G = £!4i/4i, To satisfy (6.2) for 
X < G, it is only necessary to take c sufficiently small, but the choice 
ae now depend on A, and A,. 

Hence, of the intervals (0,H), (H,2H),... contained in (0, 7’), at 
least [}477'/H] must contain points of HZ. Now, if (nH,(n+1)H) 
contains a point ¢ of Z, there must be a zero of ((4-++7w) in (t,t+H), 
and so in (nH, (n-+2)H). Allowing for the fact that each zero might 
be counted twice in this way, there must be at least 


42 
sar] > > Ale miog 7} 
9 








4H |~ 32A2 


zeros in (0, 7'). This proves the theorem. 

















TWO LATTICE-POINT PROBLEMS 
By A. OPPENHEIM (Singapore) 
[Received 30 April 1946] 


1. ConstpEr the class of ellipses H(f,A) which have centre O, area 
A, and, as a pair of conjugate diameters, the given straight lines 

S(x,y) = ax?+ 2hay+by*? = 0, 
where a, h, b are real and the discriminant 

4h?—4ab = D®? (D> 0). 

Without loss of generality we may take D = 1 or any convenient 
value not zero. For appropriate values of A there will be ellipses of 
the class with the property that on the boundary of any one of them 
lies a lattice point P(x,y) (and its image) and inside lies no lattice 
point save O, the origin. 

We seek to determine the lower bound A, > 0 of such admissible 
values of A. The number A, is plainly the same for all forms unimodu- 
larly equivalent to cf where c is any constant other than zero. We 
seek also to determine the upper bounds of A, for all pairs of conju- 
gate directions. For convenience I suppose that neither direction 
is rational. The analysis is easily modified to include the case of a 
rational direction. 

This problem, as it turns out, is exactly equivalent to the problem of 
determining the minimum (in the usual sense of lower bound) of the 
indefinite binary quadratic form f, a problem solved by Markoff, for 
an exposition of which reference may be made to Dickson (1), or the 
minimum of the product of two homogeneous linear forms. 

In three or more dimensions it is equivalent to the problem of 
finding the minimum of the product of three or more homogeneous 
linear forms. In each case, of course, the variables take all integral 
values, the set (0,0,...,0) excluded. Plainly, by Minkowski’s famous 
theorem, Ay < 4. 


2. For appropriate A there may be an ellipse of the class with 
the property that on its boundary are two lattice points P, Q not col- 
linear with O (and, of course, their images P’ and Q’) while inside 
is no lattice point except O. The problem is to find the lower bound 


A, of admissible numbers A (naturally A, is the same for all pairs 
3695.18 C 
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unimodularly equivalent to cf) and to find the upper bound of A, for 
all pairs of directions. 

This problem, it turns out, is equivalent to the following minimum 
problem for indefinite binary quadratic forms: 


Let AX?+2HXY+ BY? = F~ ef. 
To determine min |(A+ B)/2H| = M(C;) 
as F runs through all forms equivalent to cf. 
I prove that 
and that, if equality occurs, then either 
of ~ x°+-2xy—2y? = fo 
or of ~ x?+-4ry—3y? = fy. 
It then appears that 
A, = 7/(1— 2)! 
so that a< A, < 2n/N3, 
both bounds being attained. Whether A, can take any value between 
these bounds I do not know. It is plausible that there should be a 
number d < 27/3 such that, if cf is not equivalent to fy or f,, then 
A, <d. But this I cannot prove. 
It is worthy of remark that, in this problem on indefinite binary 
forms, the form x?+ay—y? has relinquished its usual dominant 


position. Naturally there are corresponding problems in three or 


more dimensions. But they are very difficult indeed. 


3. The first problem. It is a simple question in elementary 
mathematics to show that the minimum ellipse, with centre 0, 
through a given point P(x,y) with OK, OL as conjugate diameters 
(P not on OK, OL) has area equal to 27 times the area of the 
parallelogram formed by OK, OL and the lines through P parallel 
to OK, OL. Hence, if the equation to OK, OL is 


S(x,y) = ax*+ 2hay+by* = 0, 
the minimum area is W = 2/f(x, y)|/D. 
If then we consider only visible lattice points P(x, y) such that 
lf(w.y)| < 3D (3) 
(and infinitely many such lattice points exist), this minimum ellipse 


has an area less than 7. This ellipse cannot have on its boundary any 
lattice point save P and its image P’, nor any lattice point inside 
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save O. For, if the contrary were the case, it would contain a 
parallelogram of area > 2 and hence its own area > z. 
Hence the admissible values of A < z are those values of 


2m|f(x,y)|/D <2 
for co-prime integers x,y. Thus 
A, = minA = 27min | f(z, y)|/D. 
Now Markoff showed that 
lim min | f(x, y)| = 4 
and determined all the classes of forms, the Markoff forms, for which 
min | f(x,y)| > 3- 
We have therefore A, = minA < 2z/N5, 
equality occurring when 
of ~ x*+ay—y’*. 
If this class of pairs of directions is excluded, then 
Ay < 2/42, 


eq wal ity occurring when 
21 Dangz gy 
Cf ~~ x*-+ 2ry—y"*, 


and so on; the precise upper limit of Ay being 
lim Ay = 37. 
A complete statement of Markoff’s theorems may be found in 
Dickson (1). 
4. The second problem. 
[ now prove 


THEOREM 1. Suppose that 


cf ~ AX?+2HXY+ BY? 
with A= |(A+ B)/2H| < }. 
Then an admissible value of A is 
A = a/(1—A?)}. 


Conversely any admissible value of A is of this form. 

Suppose that an ellipse, with area A and centre O, and having f = 0 
as conjugate diameters, passes through the two lattice points P(p,, pg), 
Q(9;,92) and has no lattice point inside except O. Plainly P and Q 
must be visible from O. Also no lattice point can be inside the 
triangle POQ or between P and Q. Hence 


AFOQ = 3, P1%—P2h = +1. 
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Thus, by a unimodular transformation, P and Q become (1, 0) and 
(0,1), and the ellipse consequently becomes 


X?+42BXY+Y? = 1, [B| < 1, (1) 
of ~ AX*4+ 2HXY+ BY?, (2) 
A = a/(1—f%). (3) 
The harmonic condition gives 
A+B = 2HB. (4) 


Since no lattice point other than O lies inside the ellipse, it follows 
that 
|28| < 1, (5) 
this being necessary and sufficient. Both parts of the theorem follow 
from (1)—(5). 
It may be noted that, if |28| = 1, the ellipse will have three lattice 
points and their images on its boundary. 
THEOREM 2. We have 
A, = /(1— mn), 
where M = minA = min |(A+ B)/2H |. 
This follows immediately from Theorem 1, once we have proved 
THEOREM 3. We have M < }, 
equality occurring when 
of ~ «+ 2xy—2y’ = fo 
or of ~ z?+ 4ry—3y? = f, 
and in no other case. 


Corresponding to an indefinite binary quadratic form f without 
rational factors there is a sequence of positive integers 


(9;) (¢=...,—1,0, 1, 2,...). 
Write F, = [91 Gis, H; = [9-1 9i-2---; 
these being simple continued fractions, 
A,=FRH+1), B= (FB H—V/(R H+), 
Ave = HF, H,+1), 


so that eo; 
Bi+4A4; Aji, ag 1, B+ Bis, al 29; Ais, 


(x,y) = (—1)'4,2?+ B,ay+(—1)**4;,,y”. 
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The form f/D is equivalent to any member of this chain of equiva- 
lent reduced forms. For an exposition see Dickson (2). 
Consider now any member of this chain, for example 


D, = Apx*?+ Byry—A,y’, 
and with ®, the two forms 
¢’ = A, X*+ (24,+ B,)XY+(A,+ B,—A,)¥?, 
Pa (A,— By—A,)X?+(B,+2A4,)XY—A, ¥? 
obtained by the respective parallel transformations 
z=X4+Y, y=Y!; a=X, y= —X+4+Y. 
For these three forms the corresponding values of |A+ B|/|2H| are 
Ao = |Ao—Ay|/Bo, b= (249+ By—A,)/(249+ By), 
v = (24,+ By—Ap)/(24,+ By) 
since 
2A,+ By—A, = [(A—1)(Ay+2)4+1)/(H A+) > 0 
because F, > 1, Hy > 1, and similarly 2A,+ B,—A, > 0. 
Suppose that Ayo > 3, A, >A, 
so that 
2A, > 24+ By, p= 1—A,/(2A,+ By) < f, 
and it may be noted that v > }. 
Likewise, if Ay > 4, A, < Ao, then 
v4, p>. 
Hence one of the three numbers Ap, yu, v is less than $ except possibly 
when doe J, 2|4,—A | = Bo. 
Thus we have proved the existence of forms of the type required 
in Theorem 1, and also that M < }. 
5. I now show that, if the equation 
A; = |Ajuy—Ay|/B; = 3 (1, ¢) 


holds for five consecutive values of i, it holds for all i, and the only 
possible forms for which this can happen are f,, f, of Theorem 3. 


It will follow that, equivalence to f, or f, being excluded, there will 
be an infinity of admissible ellipses with A < 27/V3. There is no loss 
of generality in supposing that (1,7) holds for 

t= 0,1, 2,3, 4. (2) 
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From (1,7) by squaring and using 
Bi+44; Aj, = 1 
we get A?—A,;A;,,+A?,, = }. 
Take (4,7)—(4,i+1) (¢ = 0,...,4). Then either 
A; = Aijss (5, 2) 
or Aji, = Aj Aijse (6, ¢) 
If A,,, lies between A; and A;,,, then (6,7) cannot hold: therefore 
(5,7) holds. From (1,7) and (1, +1) and (5,7) 
0 = 2|A;,2—A;| = 2|Aj42—Aja|+2|4i—Ail 
= Bit Bi 
but B+ Bia. = 29; Atay (7,2) 
which is positive. Hence A,,, exceeds both its neighbours or is less than 
both its neighbours. Also, if A;,, < Aj, Ajay < Aj+e, then necessarily 
(5,7) holds, that is A; = Aj». 
(I) Suppose first that A, > Ay, A, > Ay. From (1,90), (1,1), (7,0) 
we obtain 
44,—2A,—2A, = 29,4,, (29g —4)A,+24,+24, = 0 
so that 
J = 1, A, = Ag+Ap, By = 24g, 
Now A, < A, so that also A, < A; and therefore 
A, = As, B, = B,= 9,3. 
Next A, > A, and so A, > A,. Therefore, as above, 
Jo = I, A, = A,+A,, 
so that, by (8), (9), Ag= Ag, B, = B,. 
Now A, < A, and therefore A, < A;. Hence 
A; = 4; = A,, B, = By = By = gy 
Since now ©, is identical with ®,, the chain has the period 4. 
Also 9,4, = B, = B, = 24,, 
9,4, = B, = 24), 


whence g,g3 = 4, so that 


9,=93=2 or gy =1, 93=4 Or 9, = 4, g, = 1, 
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giving essentially just two sequences 
(1,2) and (1,1,1,4) 
and the periods of equivalent reduced forms 
x*+- 2ry—2y?, — 2a?+-2ry+y? (13) 
[1,4,—3], [—3,2,2], [2,2,—3], [—3,4,1] (14) 


respectively. 


(II) The case A, < Ay, A, < A, can be discussed likewise. It 
leads to the sequences 


(2,1) and (4,1,1,1), 


the new forms being the negatives essentially of (13) and (14). 
[t remains to show that, for the forms in (13), (14), 


min |(A+B)/2H| = 3. 
Suppose that 
2|A+B| > |2H|, 4H?-4AB=D*?, D>O. (15) 
Then (15) is equivalent to 
A?+ AB-+ B* > }D*. (16) 
Now (16) holds automatically if AB > 0 and one of |A|, |B] is at 
least 3D. It also holds if AB < 0, provided that both |A|, |B] are 
at least 4D. Thus it remains to discuss (16) when one of |A|, | B| is 
at most 4D. Since such a form is equivalent by a parallel transforma- 
‘tion to a reduced form, we need only consider a reduced form 


AX?+ 2HXY+ BY?, |A| < 4D, 
and those derived from it by the parallel transformation 
X=X,4+%, Y=¥, 

where & is any integer. 

For (13) we then get 

A =i, 2H = 2k+2, B= k?+2k—2, 
A+B = (k+1)—2, 
but plainly, for any integer m, 
|(m?—2)/2m| > }. 
As for (14), there are the two forms 


x?+ 4ay— 3y?, 222+ 2ary—3y* 











24 TWO LATTICE-POINT PROBLEMS 
to be considered. For the first of these forms 
A = 1, 2H = 2k-+-4, B= k?+4k—3, 
A+B = (k+2)?—6, 
(A+ B)/2H| = |(m?—6)/2m| >4 (m= k-+2), 
and for the second form 
a = 3, 2H = 4k+2, B= 2k?+2k—3, 
A+B = }(2k+1)?-4, 
(A+ B)/2H| = |(m?—3)/4m| >4 (m = 2k+1). 
We have therefore proved the statement at the beginning of the 
section and have also shown that 


M(fo) = }, Mf) ms b 
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THE APPROXIMATE SOLUTION OF EQUATIONS IN 
INFINITELY MANY UNKNOWNS 


By R. C.. T. SMITH (Ozford) 
[Received 20 May 1946] 


1. Introduction 

In many of the methods used for the approximate solution of 
integral equations and differential equations the original problem is 
reduced to that of finding the solution of a system of linear equations 
in infinitely many unknowns. This system is treated by the method 
of segments (réduites), i.e. only the first m equations are considered, 
and all save the first m unknowns are neglected in these equations. 
It is this step which makes the results approximate. 

Very little is known about the accuracy of such approximations. 
This paper tries to remedy partially this lack by some work on the 
method of segments as applied to eigen-value problems. It is based 
on the standard (Hilbert) theory of quadratic forms in infinitely 
many variables (1, 2, 3). 

Systems of equations of the form 


Xa—-AY kagtg=0 («= 1,2,...) (1) 
p=1 


will be considered. It will be assumed that the matrix [k,,] corre- 
sponds to a completely continuous bilinear form, ie. one which 
satisfies equation (5). In this case equations (1) have properties 
entirely analogous to those of an algebraic system of m equations 
in m unknowns.* In general the only solution of (1) is the null solution 
in which x, = 0 for each «; but for a certain enumerable set of values 
of A, the eigen-values, there are non-trivial solutions, the x, then 
constituting eigen-vectors. 

The concern of this paper is to investigate how closely the eigen- 

values Ay”,..., A@” for the curtailed equations 


L—A s kgtg=90 (a = 1l,...,m) (2) 
p=1 


approximate to the eigen-values 4,,...,A,, of equations (1). To 
simplify the discussion it will be assumed that the matrix [kg] is 
Hermitian (kg = kg.), the eigen-values A, then being real.{ 


* (2) 94. + (4) 100. 
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The method adopted in this paper is to solve the equations (1) 
with a = m+1, m+2,... for 241, mio. in terms of 2,..., Lp. 
When 2,,,..1, Um+9)--- are eliminated from equations (1) witha = 1,...,m 
one obtains 


24— Sraptp—AS kagxzg=0 (a =1,...,m), (3) 
p=1 p=1 


where the r,, are rather complicated functions of the k,g and A. 
Obviously equations (3) can only have non-trivial solutions for the 
eigen-values A, of equations (1). By choosing m sufficiently large, 
convergent expressions for the r,g are obtained for any given A 
(because of the complete continuity of [4,,]). Bounds are readily 
obtained for the r,,, and the maximum difference possible between 
the eigen-values of the exact system (3) and those of the approximate 
curtailed system (2) can be estimated. To show the convergence of 
these differences to zero, with increasing m, a slightly stronger 
assumption than complete continuity seems necessary. In addition 


it is assumed that 
m <<] 


lk ,|2 
> [kag 
a=1 P=m+1 
tends to zero with increasing m. This condition is satisfied, for 


example, if 


converges. Since 


where the two sides converge or diverge together, the convergence 


a 
continuity, which implies merely that 1/A, tends to zero with in- 


creasing y.* 
Further slight extensions of the results to systems of the form 


Ly+ > haptg—A > kuptg=90 (a =1,2,...) (4) 
= f= 


of > |k,g\? is a slightly stronger restriction than that of complete 
,B=1 


will be indicated and an example treated. 
Throughout it will be assumed that (x, 2,,...) is a point in Hilbert 


«oo 
space, i.e. that > |x,|* converges. 
a=1 


* (3) 1558. 
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2. Preliminary results 

This section will give the most important theorems needed later. 
Proofs are to be found in (1, 2,3). Throughout, the parameter A will 
be assumed to be real. 

The definition of a completely continuous bilinear form > k,g%,%g 
is that it satisfies 


| > kapr.¥g— > kapta¥s| <€, (5) 
a,B=1 a,B=1 


for arbitrary positive «, for m sufficiently large, and uniformly for all 


Ty, Yg satisfying * @ 
2 Ia? = 2 |vel? =1. 


Then s kaprtaYp| < €, 


a,B=m+1 
where > Itel®?= D> lygl?=1 
a=m+1 B=m+1 

since the left-hand side is a possible value of the left-hand side of (5). 
The maximum value of 

= « 

lepaasatae™e¥e| 
under the conditions 

> |z,|* = t, > lyg!* =1 
a=m+1 B=m+1 
will be known as 6,,. Therefore, for completely continuous forms, 
bn <€ (6) 


for arbitrary « if m is sufficiently great. 
{t is well known* that equations (1) can be reduced to the form 
Yy—(A/A,)y, = 9 (y = 1,2...) (7) 
(where the A, occur in order of increasing magnitude) by a unitary 
transformation 
Yy = 2, Way Tq (8) 


1 Y 


The defining property of a unitary transformation is 


= a 0 ( al B), 
aor by = Sag» where Sap = ( (x = ), (9) 


* (3) 1562. 
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which implies 
By using (9) the transformation (8) can be inverted, giving 
z,. = s Wey Yy (11) 
y=1 


It is obvious from (7) and (11) that 
(12) 


Loy = Way Vy 
is the eigen-vector corresponding to the eigen-value A,. Comparison 
of (1) and (7) gives 

ww 
beg = > Corey, ass 
y=1 Y 
which is equivalent to 
=. F< D ys, Week (14) 
, bead of ay “Be aB> 
Some results on the solution of non-homogeneous systems of 
equations by iteration will now be derived. They are, of course, well 
known. Consider the system 


ta—A> kag =Cy («= 1,2,...). (15) 
p=1 


By the unitary transformation (8) this is reduced to 


y(1—A)Ay) =, (y = 1,2), (16) 
where 


Equations (16) have the solution 


, 
Cc 


y 
/ 


oe Ta,” 


and thus equations (15) have the solution 
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directly in terms of the kg; but, if |A/A,| < 1 for each y, 


20 ®py+ 2 "ath 


Continuing this process gives 


Way WBy __ s ae 
p (A,)” Nyse Inn =1 an, 0 12"** Mer B 


== = k&, say. 
It will be seen that 


S) War By — = Sgt > rey. 
n= 


—NA, 


_ Thus (17) can be cevniil in the form 


_ S \npn) 
a = 2, clot 2 | 


(20) 


provided that |A/A,| < 1 for each y. Since the |A,| are non-decreasing 
functions of y, this is equivalent to |A/A,| < 1. The method used in 
deriving (20) requires the assumption that [kg] is Hermitian. How- 
ever, the result is true for arbitrary kg provided that |A|b < 1, where 
b is the (Hilbert) bound of [/,g].* Most of the theory in this paper 
follows from (20) and could be extended to k,g with no symmetry 


properties whatsoever. 


A few further remarks will be made. The first is that the (Hilbert) 


* (2) 92. 
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bound of the matrix [k,g] is |1/A,|. This follows sincé the bound b 
is the maximum value of 
> kp XOX), 
a,B=1 
where the X®, X® are arbitrary apart from the conditions 


@ is 3) ( 
2 xe? = P» |X| i t. 


« Pad 
But kgpXOX® — 
wha 1 - 


where the Y’s are related to the X’s re the transformation (8). It is 
easily seen that 
0 oe a 
YXOR= SOP, FS |XPPe= > Vo! 
a=1 y=1 B=1 yt 
from the properties of unitary transformations, and it is then obvious 
that the maximum value of 


is attained when 
\Y@| = 1, \Y@) = 1, 


YO=0, Y2=0 (y>1) 
and is equal to |1/A,|.. Thus 
b = max Ss kop Xwxp| == |1/A,|. 
a,B=1 
By a similar argument upper and lower bounds are obtained for 
the value of the Hermitian form 


~ | . npn) | & 4 
Pes [Suet 2.4 kp |*aXp, 
namely 


. n (n) | ¥ x a 
(Bug+ > kp \| Xa BS TP 


(Sap-+ 5 rk) XX < 
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This result is obvious since 
a) 5 Cs r ; Ax bad VY. 
+ > Avkey) XX M21 aah ‘ 
Pa ad p> af | —- 


[f all the eigen-values A, are known to have the same sign as A, the 
coefficient 1/(1-+-|A|b) in the first term of (22) can be replaced by 1 


The bound of the Hermitian matrix [5.p+ s dnd] is 
n=1 
1 1 
= 23 
1—|A/A,|  1—|A|o° ei 
This can be used to examine the order of magnitude of the x, given 
by (20). For 





Z_ = C+ “a Ce >. Anke 


2) 


A> bolbat SHB) ep 


ee 


Ais {3 Neel} 
so that is,| < [e,|4-—— =a —, (24) 


3. Reduction of the infinite system of equations to a finite 
system 


The system of equations 
t_—d¥ kag =O (a =1,2,...) (1) 
p=1 
is readily reduced to 
m m 
ty— ¥ Tepxg—A > kagxzg=O (a = 1l,...,m) (3) 
p=1 p=1 


by solving for 2,41, %m4+9)--- in terms Of 2},..., Z,, by iteration and then 
eliminating. To do this, equations (1) are separated into two groups 


Ly—A s kaptp—A s kygxp = (x = 1,...,™), (25) 
p=1 B=m+1 

Ly—A s kup&p =A 3 kup&p (x =: m+1,m-+2,...). (26) 
p=m+1 p=1 


To avoid confusion aad last group will be written as 


Xy —A s Kim « = A 2 Papp (ax = m-+-1,m-+-2.,...), (27) 


B=m+ +i 
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ities ka = : a 2008 m), 
-s CoB (a, m) ° 


Equations (27) can be solved by iteration, giving 


ta 2%, (Port BK) S bypag (a= mtLm+3,.), (28) 


y=m+1 
the bound of [K,,], is less than |1/A|, i.e. if 
[Alby <1. (29) 


provided that b 


m? 


From (6) this is satisfied if m is sufficiently large. On intuitive grounds 
it might be expected that b,, is an approximation to |1/A,,,,| and 
therefore that |A\b,, < 1 iy |A| < |A,,|. It will be assumed that 
(29) holds. 

If the unknowns 2,41, Lm+2)--- given by (28) are substituted in (25), 
there results 

a.—\S kom—¥ SF &[0.4 See Sk O 

e p> of B il ay | wt 2 ¥ 2, B B 
(a = 1,...,m). (30) 


This is of the form (3) with 


@ ie) 
_ l z(n)\ J. Je 
rep =O eat See | ay hep aia 
Obviously r,¢ is Hermitian for real A. Equations (30) together with 
(28) are completely equivalent to equations (1), provided that the 
convergence condition (29) is satisfied. 
Before passing to the next section it should be observed that the 
result of applying (24) to (28) is 
; oo — co 
al{ > Saal { > 
‘y=mt1 y=mtl 


—|A\b 


s 


kyp | } ry 





> =1 





la| < \A Ss k 
p=1 


m 


‘ ‘ Kayl*|'{ > > Il} 
< ail ial?) saa Nts ial 


(a = m+1,m+2,...), (32) 
<{X laal*}"{E lbal?}" has 





where Schwarz’s inequality | 
a 


been used in simplifying the result. 





: 
| 
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4. Comparison of the curtailed and exact systems of equa- 
tions 


For the curtailed system of equations 
ty—A ¥ kagxg =O (a= 1.,...,m) (2) 
_ A 


in all the results given by equations (7)-(14) summations from 1 to m 
replace summations from 1 to 0. To indicate that the curtailed 
system is being considered the symbols will be modified by the 
addition of a superscript (m). The m eigen-values A{”,..., A% of (2) 
are found from the equation : 


det (5,3—Akag)™ = 0, (33) 
which is the condition that (2) has a solution other than the null 


solution x, = 0 (a = 1,...,m). If wf) is the normalized eigen-vector 


corresponding to A, equations (2) can be reduced to the form 
Yy—(NAM)y, =O (y = 1,...5m) (34) 


by the unitary transformation 


™ 
Yy = 2 Oey a 
- (35) 


m 
(m) 
2, > wey Vy 
y=1 


For the exact system (3) the results are not as simple. The 
determinantal equation 


det (84g—Akag—T ag) == () (36) 


is the condition that equations (3) and hence equations (1) possess 
a non-trivial solution. It will be shown that (36) determines m eigen- 
values A,,...,A,, to which A{™,..., A are approximations provided that 
the r,g are small. The values of the x, (« = 1,...,m) will be propor- 
tional to W,,..., Wmy, the first m components of the eigen-vector w,, 
of (1). The remaining components are given by (28). In general the 
vectors {Wyy,-++5 Wy: 0,005 0,...} will not be orthogonal. 

[t is convenient to transform equations (3) by (35) rather than to 
discuss them as they stand. This gives 


Ya(1—A/AM™)— ¥ rigyg = 0 (a =1,...,m), (37) 
B=1 


3695.18 D 
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where 


m 
Be is > (I) » ,,(m) 
7 = Ww Ww r. 
a 
4 “_. 


eo oO m 
= r2 > (S,0+ > Ky > kya kine Om we 
n=1 7 ye=l1 


7,0=m+1 


> (S0+ > Ky) rhe) (38) 


7,0=m+1 


™ 
and rio) — > kota (39) 
= 


Like [r,g], [vag] is Hermitian. By applying the unitary transforma- 
tion which reduces [K,g] to diagonal form (compare the derivations 
of equations (21), (22)), 74g can be expressed as 


% ‘ c ‘ 
e- 2 oe- 2 poe 

yt t—M tly, AX, Aly)! (Alin) 
where fini, Mm+g-+ are the eigen-values corresponding to [K,g], 
|A/ey| <1 (y = m+1, m+2....) since by, = |1/Hmii|, and 
<7 | rio |2, 

d=m+1 

It is easy to see that ITepl® < Tee MBB» (40) 
with the equality sign holding if and only if x) = kac\®) for each y, 
where k is independent of y. It is also obvious that rj), is a con- 
tinuous function, zero for A = 0, and increasing as |A| increases, since 
s contains a factor A. 
Applying the result (22) to (38) gives 


@ 


ao 
| pla) |2 (cx) | 2 
Yr > |r| 
y=m+1 <x” < yY=mti1 


i+, = = 1—Al, 


m 
A rough upper bound for r,, is given by 
m 


7,0=m+1 
RS > [54+ FrmKol kok 
—_ 22 + nK n \k k i 
y=1 steak \ W Py 7 7 yn 
For 7’ we have 


@ 


r2 > py [qyn |? r2 > > [yy |? 


y=1 y=m+1 “yr < y=1 n=m+1 


+], Jl, 
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™m oO 

The assumption that 5 >  |k,g|? tends to zero with increasing m 
a=1 B=m+1 

ensures that r’, and hence r’,,, also tend to zero. If the conditions 

of the following theorem are satisfied for all m sufficiently great, it 


is necessary that r/,, > 0 in order that A%” + A, as m > 00. 
THEOREM. The eigen-value d, (a = 1,...,m) of the system (1) will lie 
in the range 


1—(¢_ nse dk SS 
Fondant + 1—P'/taJ yen 


, t ) 
A/AM™ < 1—(Th4),= Fin r |e 43 
< al i < ( aa) A=) 1 1" ta} janom ( ) 
where X', X” are any numbers such that 


4 r’ /t ’ : 
] (rpg) yay 1 Pe < A [ag < 1—(ha)r-r{] + 
| i—r [tad n=nom 


1" Ite 
1—r’/t,J’ 


r'|t 
=. . <X"/A™ <1, 
( xa A=X | 1 r" |t — SS 1a 


and t, is any constant such that 
0<t, < min |1—A/Ay” | (B A a; 1—2Ar.a)aanom < A/AM™ < 1) 
provided that 
(r’/t)a=am = 4. 
If X, is the smallest positive eigen-value or the smallest in magnitude 
of the negative eigen-values, so that 
1—A/AQ” (BA a; 1—2Area)a-agy < A/A” < 1) 
is necessarily positive, (43) can be replaced by the smaller range 


vy} 
1—(rea)aer} rans < AglAM™ <1 —(rya)pey (44) 


1—r'/t, A=Am™ 
where X', X" satisfy the slightly modified conditions 


‘ r’ /t 
1—(Tya)a=rvm™ 1 +i-vh } xm 
l"aJA=XZ 
1" Ite 


< *1)(m) a tas 4 ae Ay ’ 
= A [AZ —_ (Taa)a=r | TI—rht, A=Aw 


1—(readrer SAA <1, 
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again provided that (r’/ta)a-am < F. 
This theorem will now be proved taking, for definiteness, « = 1. 
Equations (37) will be separated into 
m 


yy(1—A/Ay” —141) a rip YB > = 


p =2 


m 


and Yoa(1—A/AL-)— ¥ rhigyp = ay. («> 
p=2 


Dividing (46) by (1—A/A®”) gives 


Tap¥p_ _ _tath_ (4 > 4) 
2 1—A/A) 1 — AJA” 


and this system can be solved for y, by iteration. Substituting the 
result in (45) reduces this to the form 


Yy(1—A/AY” —14, 1) = 9, (48) 


nt ra 


— 9% — TT nugTR 
where e la! ta ap” pi +-.... 
ia Duca A NG +20 A, Ax” )(1—A/Ag”) 


It is easy to verify that e,, is continuous and real for real A. The 
condition that (47) can be solved for y, by iteration is that 


ror 
a ae eo 
[re 

2\4 


- * yr m | 
Ow | ‘a | > |S A/Aw) | j ° 


| 
Let ¢ be a positive constant such that 


-e l \ rm | (a a ae I(r" m r Am) < 1 


\. es frotie| 


If r’—r}, < t, the formal work used in deriving (48) is certainly valid. 
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Bounds will now be sought for e,,. We have 


, al Ig fu 
a bare aes say +20 —A/A™| 1— Apap 


—r 
—- a , where 1 —2(r}1),-yom <A/AY” <1 
m) 
adie (49) 
since r’—r}, increases as \/A{” increases. 
If (1 —Fa)rengo < HF, 
then lex] <1, Where 1—2(r},)) nm SA AM < 1 
It will be assumed that this condition holds. 
If A, A’, A” are such that 
1 —2(r5a) ranger < NAL <AAM <AAM <I 


then 1—(r}3)) m +n ot) | 
| t t Jajam 


< (l—1—e11)a 


’ r—ry r’ —13\) 
i—fr = n/(i— My 
oe t t } Jaa 
Further, if 


N/A < 1—(rhy)p_yef /(: Sarna || 
| t t A=A\™ 


Q fg! y’ —r, 
ond "Am > 1—(ri3), y fy _# “u /(1 a u)) : 
| t t A=\™ 


then (l—r},—e,,—A/Af”) is not positive for A = A’ and not negative 
for A= A”. From continuity there is a zero A = A, of 
f /)(m) 
(I-13 —A/Ay””) 
in the range 


Yr’ —r r’—13\) 
1—(r3,)). an ee ee 
(Tida=x | a; | r owe 


eee eS 1] . (50 
> r rj > (711)a=r | t t A=aAm™ ( 
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The effect of replacing (r’—r},)/t by r’/t is to increase the range and 
gives 
, r/t , 
]— (734) rfl+ | <A a < l- (Ty), x{] 
7 | r a” 
where A’, A” are any numbers such that 
rit \ 


2(r; wy < AA < 1—(75,))~ 90 ly + ——__. 
1)d AY ‘Tis | ( Wa nN | 1 1—1"/tJ yy’ 


1—(r)aex(I—-—| <A <1. 
(Ti)a *) r/tJ aman =~ 1 


By taking A” = Aj”, it is obvious that the condition 
x’ It | veo 
=F < A’/A™ 
A AY” 
can be added. The first part of the theorem has thus been proved. 
The second part is deduced in an exactly similar manner using instead 


— (Ti)a=a 1 - 1. 


of (49) the inequalities 


z ne r’- J y’- J e 
0< an < eesti a 


which hold if (1—A/A%”) is positive for each a. To prove (51) it has 
to be shown that e,, is positive. Now 

m. Ine? 12 m. J J J 

€ = S a! + ~*~ Tia" oap "BI Bicvg 

av 1—A/AM™ (1—A/AM%™)(1—A/AY™) 

x=2 ” x, B=2 ibis P 
and the sum of the terms on the right (excluding the first) is less 
in magnitude than 
Mm , , 1 , 
. ite Tea)" S " gp)" rp fy. Py = 

| A Mm Ly =i Ag” \ <a Am ' 

aX g=% y=2 Y 


2 p=2 


m 
Z ry x 
ao aA aS” 
x 


, ° , , 
Nal” (r’—r},)/t 

ANA” 1—(r’—714,)/t 
| i 2 
"la 


1— AJA) 
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Since the first term in the expression for ¢e,, is positive and the 
remaining terms have a sum less in magnitude than the first term, 
€,, is certainly positive. Having established (51) it is then simple 
to derive (44), ie. the second part of the theorem. 

If the bounds are to be of practical value r’/t must be small, 
preferably much less than }. If 7’/t is not small, a more careful 
discussion is needed. This case corresponds to two or more of the 
A” being nearly coincident. Suppose that A{” ~ AY”, so that t is 
small. Then ys,..., ¥,, are eliminated as before, yielding the equations 

yy(1—A/AY” —1731— C11) —Yol(Ti2 +12) = 9 | (52) 
= Ya (21 +21) +Y2(1—A/AY” —T92—€ 92) = 0, 7 

ry r’ —1i—To9 
1—(r’—1}, —Toe)/E ‘ 


where 


a 
"22 
, , 
1 —(r’—1r,,—1o9)/t 
’ 1 m L 
(731)° (729) 
1—(r’ —14,—19)/t 
and ¢ has the modified definition 
: (m) . ie He’ te’ 4 (m) \(m) — 

t 1—A/AM™ | (a > 2; 1—2(rq1+192)y=yem < A/AY™, AJAY” < 1). 
Equations (52) will not be discussed in detail. However, neglecting 
C11 y9> 01> €gg and using the result |7r}.|? < 7}, 79, it is easy to show 
that there are zeros A,, A, of 

(m) od f 
| 1—A/Ay”—1y —Tis 
| , 
—Tm 
in the ranges —1y > A /A™ > 1—1yy—12: 


and 1 >A, /AM > 1—T 9. 


(53) 


5. Extension of the theory 
There is no particular difficulty in extending the theory to systems 


of the form 


Lat Ss haptp—A > kygxtg= 90 (a = 1,2,...). (54) 
p=1 p=1 


The same conditions will be placed on the h,g as on the k,g, namely 
that [hg] is Hermitian, completely continuous, and that 


m oO 


1 ag” 


a=1 B=m+ 
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tends to zero with increasing m. The further assumption that 
[5.g+,g| is positive definite will be added. 

By a unitary transformation 


wi 


the equations (54) can be reduced to the form 


ox 


YA1+h,)—A > kigys 
j= 


Writing y(1+h,)' = yi); 
brings (56) to the form 


Yo—-A > kegyg=0 (a= 1,2....). 
1 


p= 


Since [5,,+-/,g]| is positive definite the (1+-/,) are positive and kj, 
is Hermitian. The eigen-values of (58) are therefore real. Finally 
the unitary transformation 


reduces (58) to z,(1—A/A,) 


The eigen-vector x,, corresponding to A, is given by 


x 
= bs iy UnB WBy ™ 61 
<: (1+-hg 2? \ ) 
Exactly the same process as that used in § 3 can be applied to (54). 
This gives 
m m 
+ > hygxg—A> kygxzg=90 (a= 12,...,m), (62) 
ee fa, Pn 


Pp 


where 


(B,<+ >" K)} (Niky Nay) (Akeg—hep) (63) 


(y,e >), 


(y,e > m), 
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and the K%? are defined as before. The properties 


7g |* San TBB> 
> | May—Iray|* > Aky— 
=m+ <r < y=m+1 
iF) AlDn a 1—|A\b,, 


still hold. Minkowski’s inequality then gives 
| S ip tal ¥§ i2\*\? 
{A( >  |hayl?)'+( > [ray l*)'| 
y=m+1 y=m+1 : (64) 
1—|A\b,, 
Applying the transformations for the curtailed system analogous 
to (55), (57), (59) reduces (62) to 
m ‘ 
ie | ] —x, Am”) - Tapp == @ 
=1 


=(m) mm) (m),,(m) 
Urn U P Ug Ugg 


where ro = - r, ; 
ap ( +H) (1A)! 


y-€.7,0 =] 


By using (61), 7g can be rewritten as 


m pm) em) 
, - © 
2 = r, = : 67 
a8 2 ” (P= )(2), ( ) 
Ry 28 


y-€=1 
and finally it can be brought to the same form as (38) (which it in 
fact includes as a special case), 


x 
e) roy =(B) 
op = > 80+ > A" KY ie 
7,0=m+1 
oy 


Zim)" 


where i = s Ake —hy 9) — 
y= 


is , on 
For 7, the inequalities 


. ( ~ ( 
| p(x) | 2 (ax) | 2 
> r9 | > 1 | 


d=m-+1 a 6=m+1 (70) 


i+fe, “= 1—(AM, 

hold. 
The rg have the same properties as in the previous section and 
exactly as before the eigen-value A, of (65) is determined by (43) 


with o 


f=m> ft. 
a=1 
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Bounds for r’ are given by 
ati 
x 
(1 +h min 
6. Practical application 


Usually in (43) 7... will represent a small correction to the eigen- 
. Xa 2 


value, the terms 


on the left and right can be neglected, and it is sufficiently accurate 
to evaluate 7, for A = AY”. The problem of estimating the propor- 
tional error in the approximation A‘” to A, is simply to find bounds 
for 7’,,. To be definite the first eigen-value A, will be considered, and 
to simplify the notation the superscript (m) will be dropped except 
from Aj”, as the Z,, %yy, Wag, Wag, Aq all refer to the curtailed system. 
It will be assumed that the eigen-vector x, = x,, has been calculated 
but not the w,g and wyg. 

To use (70) the value of |z,| is needed. The best result to be got 
from (61) without using the values of wu, and w,g follows from 


p 


x 


m ee w 2 
7. r - = Be Bi 
heat 


— l+hg 


1 


x 


m 


It is (1+hez) 


B/min & 


This can be used to rewrite (70) as 


m 
where Ro = > (Myp—h,p)&,. 
y=1 


/ 
If |h,.g| is positive definite, the hg are positive, and 
ut 


(2g) max : 2 hs > higp. 


Q 1 


5 
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In this case (71) can be replaced by 
m 1 m j \ + 
|? < |z,|/? < (1+ bee) > |x,/2, 74 
2, | —— ( + 8) 2, ‘ (74) 
m 
and (hg)max: (hg)min in (72) replaced by > hgg and 0, respectively. 
B=1 
From (73) 
dS Fx Y Wkwkgo—Akawhgo—Ah okgothaghge). 


a x K B 
m+1 
(75) 


R,|? = 
1 


a,p==1 
[t is convenient to consider separately the various terms on the right 


of (75). For the first term 
s .2ltf li: 2\? 
ho “BA 

6=m-+1 S pa B 


=m+1 


(76) 


m 
|_)2 x_| aol! 
xF p- ” Silline Bi 
A similar result holds for the term 
$ z,2, > hgh 
a2 “aX b,6 /’B@- 
0,8 1 . By m+1 Fr 
Often the result (76) can be improved. Thus, if /, vanishes for 
x — 
@ odd (even), then , > k,ekge vanishes for «+8 odd (even) and 
m 
the terms with a+ odd (even) can be omitted in the upper and lower 
For the term 
n a 


-r S XX > ko hge 


X 


bounds given by (76). 


i ' 6=m+1 


Schwarz’s inequality gives simply 
m J 
-A 5 Ta™B > 


x,B=1 m-+ 
. ty oS 2\! 
ko | > hye! je (77) 


6=m+1 


Again (77) can often be improved. If k,9 vanishes for «+6 odd (even) 
and hgg vanishes for 8+-6 odd (even), the terms with a+ odd (even) 
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can be omitted. Similarly, if 4,9 vanishes for «+6 odd (even) and hgg 

for 8+-@ even (odd), the terms with «a+ even (odd) can be omitted. 
Use of (76) and (77) in (72) enables bounds to be found for 7}. 

If there are no simplifying circumstances, the upper bound 


m 


(2... [ia f wou elty al (78) 
a / t 1 b ic 
Fe? ” \g . 1 : 


is obtained for 


x 1 1 9 
. 12\! . 2\2]? 
> koe | T , a , h x0 | ° 


m+1 

Similarly a lower bound is 
m x 
ae 3 |x? > 
a 


x=1 6 1 


~{¥ |e > lhl? +{ ¥  |hggl?}*]. (79) 
| 2 | 2 2d lhael*} |p (F 
x=] 1 
For calculations the values of 
| \! ( ¥ \! 
i > |hal?} 
é “ad ? : x0 
\o—m+1 \o—m+1 
and b,, are required. The first two will be simple to calculate either 
entirely numerically or by first obtaining an analytic expression for 


Determining b,, may be more difficult, though the result 


i 


b. : s k |)" + S 
, m 


2\! 
/ xp ' 
7 x,B=m+1 


xB 
1 


is useful. In general |A,\b,, will be small and may be neglected in 
comparison with | if it cannot be estimated. 
An alternative procedure to the use of (72) is to consider the 
determinantal equation analogous to (36), 
det (O4g—T ag th 18 -k yg) = V0. (80) 


If (80) is expanded and upper and lower limits are substituted for 
the r,g, making full use of any special properties the r,g may have, 
more accurate estimates of the error in the approximation A\” will 
be obtained. This is because (72) suffers from the uncertainty in the 
value of |z,|". However, expanding the determinantal equation (80) 
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involves more labour. In the case [hg] = 0 the two procedures will 
give almost the same results. 

Example. Leggett obtains the following system of equations* 


4 
r?+m?)?A,—imQ pa C,,A,+PR il A,+ 


(r?-+-m?)? * 


Z Es 
PRL, ¥ «, A,+U, ¥ (—1)", A,,| =0 (r=1,2,...) (81) 
n=1 n=1 

in discussing the elastic stability of a long and slightly curved strip 
under shear. This system has a solution (A,, Ag,...) with not all the 
A, zero for certain critical values of the shearing stress Q, the smallest 
of which corresponds to the buckling load. Leggett treats this system 
by the method of segments, comparing the approximations to the 
buckling load given by the second-, third-, and fourth-order 
determinants. 

In (81), m is a parameter depending on the wave-length of the 
buckle, and PR depends on the curvature. To bring (81) into the 
notation of this paper, replace m, r, n by p, a, B, and write 
(u2+o07)A, = 2x,. Dividing (81) by (u?+.a?) gives the form (4) 


x a 


+ > h, x Fk re=0 (| 
- ap B . hed) 


with hag = PRH3+ PRHY, 


2ap[1—(—1)*+F] 
= *§ = ms 
kag (a? — B?)(u?+-0?)(mP + 


a 
HY? = — 84g, 
af (?+-0?)4 ap 
H®) - [1+( a 
yn ee 


a 18 x2)3 
and a, b complicated functions of ju. 
The work necessary in estimating the error of the second approxi- 
mation (from the third-order determinant) will be outlined. This error 


* (5), 17 (52). 
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can then be compared with the difference between the second and 
third approximations. 
It will be seen that [h,g| and [k,,| are Hermitian, and that h 

af af aB 
vanishes for a+ odd, k,g for x+f8 even. By (76) and (77), using 
the simplifications due to the special properties of hog and k,,, an 
upper bound for 

> |i? (m= 3), 

d=m+1 

is the following expression with the upper signs throughout and a 
lower bound is the same expression with the lower signs throughout. 


x x 
| x, |? > |kyel?-+1 a2)? > 
6 


1 


Of course, if the lower bound given by this expression is negative, 
it should be replaced by zero. 


The calculation of 


x x. ” 

is simple. Write kyo\2 as > |kye|?— > |kyo|?. The quantities 
ti 6=1 6=1 

kg| (a, 6 are known from the work on the curtailed 


x 
system. For > |k, |? it is easy to obtain the result 


‘ oth Aap 3u2—a2 
Ba a) __ 60 h dap ~ ¥ (a odd), 
(u*-- x?)4 sinh? $a Larus peta® 


ss l tanh dap 3y?—a? 
= 2 ee ae =| (x even) 


' cosh? $a 


1 
(u?-+ a2)! br peta? 
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i 2) 
by standard methods of summing series. For > {hyg|? 
6=m+1 


x 


S lhael? = 4(PR)Ya+6(—1)297(H,)2_ > ca “Hy) 


6=m+1 6=m-+1 2 


~ 


4(PR){a—b}(H,)? ¥ (Hy)? (a odd), 


odd @>m 
4(PR){a+bP(H,)? -> (Hy)? (a even). 
\ evend>m 


It is then best to calculate 


> (A)? 


even@>m 

numerically. 

Further (hg)max> (2g)min (8 = 1,2,3) are required. These are the 
createst and least of 11, hoo, wg, Where 41, 3 are the roots of 

—p+hy hys a 
hs —p+hgs| 
Of course, it is easy to see that [h,g| is positive definite, so that the 
cruder estimate (74) can be used. For 
m 


hptatp - = PRS HQ »)?+2PR(a—b)| > Hata] + 


odd a 
m 2 
+2PR(a+)| S H,%,| 
evena<2 
and each of the expressions on the right is positive. 

Finally the value of |A\b,, is required. Let [a,g],, be the matrix 
[ag] with the first m rows and first m columns removed, and let 
bla.g|m be the Hilbert bound of [a,g]|,,. The more important 
properties of bounds are given in (1, 2, 3). We have 

|A|bn iat b[hap—Ak a hn 


< bfhaglat A|b[ hag bmn 


< PRU HB| m+ P ROH Bm + \Alblkaphm: 


Now b[a,g] is the maximum value of 


s Apr yy 
p= 


1 


a, 


where 
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If ag = 4g, and d,g is real, the maximum is attained for x, = y, 
and y, real. For [H%2|m> 0 bf H® 


aplm 3 is the maximum value of 


| p(x)? | 
| a=m+1 (u?+-0*)4| 


under the condition 


Obviously OAS Im = = 
For [AB )m> we have 


b[H2y],,—max| SY  2(a-+b)(Hyx,)(Hgxg)+ 


m 
evena,B>m 


+ > 2(a—b)(Hy2,)(Hp xg) 


~ 
odda,B>m 


0 | 
max; 2| a +6) | zz (HX ,)(Hg xg) | + 
evena,B>m 


+2\a—b|| > (H,x,)(Hp2,)}} 


odda,B>m 


< max 


. 2 
(x4) | 
odda>m 
: max x||§ a: 
evena 
For [k,g],,, more difficulty is experienced, and appeal is made to a 
series-inversion formula due to Titchmarsh (6). This states that, if 

- i= Ug 
W,=- 5S B 


ss T gw atB+h’ 
Tha ie 2 
then 


and 


Some preliminary transformations will first be made. We have 


2apf(—1)*—(—1)] 
1 2 PNG aE LE 


y 


(—1)*x, yp. 


R 
X,p 
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Writing (—1)*x, = 2,, we get 


$ aja =| 5 eee 


Yg) = - Za Yp\- 
apes OPM =| 2 Pe Ne 


The expression on the right is a bilinear form with symmetric real 
coefficients, and therefore attains its bound for z, = y, and y, real. 
Thus . 


|< ‘ —(—1)] | 
b[kagln = max} ; in a A ae” A 
is Bol Fat) (ue + BH) 


where 
’ 
Since kyg vanishes for a+ even, this can be written as 


7 2a8 | 


blk - . 
Pe (PB) +a) +) 


8|m = 4max| 
jevena>m 
odd B>m 


> amet 


evena->m 
odd B> >m 


pa eae 
evena>m (u?+-0?)(u?+-B?) r+ B 
odd B>m 


1 
+g) uae 


Let a = 2y, B = 2e+1, WU, = Yos Ve = Yep, 


ps mel 
(a? —B*)(u?+-0)(u?+-B?) 


svena>m 
odd B>m 


yily Ve 
>sm (uP 4? M ee +( 2e+1)3\y 
ies 


Ei = U,, 
p+ 4y" 


2) 
1 
Then blkgln = 4 max ZZ U, A + ne f 
[ af] Fn ye ytet+4 —y+e+4 


€-im 


Let 


E 
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Defining U, = 0 (|y| < 3m), U_, = U,, we have 


bl hap lm = 4 max’ p23 »2 " U, I 


“~e-+3 


- 
€ 
e>im y=—@ 


4 max| ¥ V. W, 


‘emtm 


< 47 max 


where 
and y W: 


€= 


The condition 
gives 


which will be written as 


Obviously 


bkaglm < las cna! ’ Ups 
0 (u*-+ 4y*)i pe" T (2 T 1)*5_ max 


| 
eile 2 _| { : 
{2+ 4/"}max \u?+(2e 


It is thus simple to evaluate |A\b,,. Substituting the lower bound 


a 
xX 


a 

for > ||? onthe left of (72),and the upper boundfor > | Rp)? 
6=m+1 d=m+1 

on the right gives bounds for r},. 

Table 1 gives upper bounds (u.b.) and lower bounds (1.b.) for 7}, and 
upper bounds for r’ for various values of PR and yp; also upper and 
lower bounds for the fractional error. The upper bound for 7’ is 
calculated from 

m 
» 


Jv x=1 
, i. 


an 


= | hag 212 kop 2] 


(1+-hg) nin ™ (1+-hg)min( 7 A\b,) 


The result (44) is applied to obtain the extreme possible values of 
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the error. The eigen-values A,, Ay, A, for the third-order determinant 
are A,, —A,, 0 and therefore 


t= min|1—A/A,| = 1 («a €1,0 <A/A, < 1). 


TABLE 1 





L.b. u.b. u.b. 1.b. u.b. 
for Tr}, for ri, for r’ for error | for error 








0-794 0-0036 0-0117 0-0300 0-0036 0-0120 
0-911 0-0023 0-0170 0-0433 0-0023 0-0178 
1-072 0-0009 0-0295 0-0718 0-0009 0-0318 
1-162 0-0007 0-0483 0-1065 0-0007 0-0541 
1-245 0-0013 0-0728 0-1481 0-0013 0-0854 
1-323 0-0024 0-1213 0-2215 0-0024 0-1558 
1-323 0-0030 0-1801 0-3068 0-0030 0-2599 








ioe) 
A direct arithmetical computation of > 2% determines closer 
6=m+1 
bounds for rj, and thus closer bounds for the error. The results are 


given in Table 2, the same upper bounds for r’ as in Table 1 being 
used to calculate the bounds for the error. The fractional difference 
between the second and third approximations (1—A{"/A{*) is also 


included. 
TABLE 2 





1.b. u.b. 1.b. u.b. 
for r, for ri, for error | for error | 1—AY/AP 


0-794 0-0037 0-0040 0-0037 0-0041 0-0023 
0-911 0-0027 0-0037 0-0027 0-0039 0-0026 
1-072 0-0022 0-0044 0-0022 0-0047 ——- 
1-162 0-0031 0-0088 0-0031 0-0099 — 
1-245 0-0051 0-0194 0-0051 0-0228 0-0106 
1-323 0-0081 0:0464 0-008 1 0-0595 — 
1-323 0-0099 0-0868 0-0099 0-1253 0-0466 








Closer bounds for the error can be obtained by changing the 
variables in (82) to 
Yq = £,(1+PRH®)! 
and dividing the ath equation by (1+ PRH®))! or, alternatively, to 
the variables 


= 


Yy = %,(1+PRH®Y)! (ax 


Ya = Ly 
and dividing each of the first m equations by 
(1+ PRH®)! (a = 1,...,m). 
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The second procedure has the advantage that most of the mani- 
pulation and the arithmetic used in obtaining Tables 1 and 2 
can be taken over unchanged. Table 3 gives the results from this 
method. Again } Rj was computed directly. An exact value is, of 
course, preferable to upper and lower bounds and is perhaps simpler 
to calculate. 


TABLE 3 





1.b. u.b. u.b. L.b. u.b. 


, 


pe Jr 7. for r; for r’ for error | for error 


-794 0-0037 0-0040 0-0300 0-0037 0-004] 
‘911 0-0028 0-0036 0-0423 0-0028 0-0039 
-O72 0-0026 0-0041 0-0670 0-0026 0-0044 
-162 0-0042 0-0079 0-0952 0-0042 0-0087 
-245 0-0083 0-0170 0-1284 0-0083 0-0195 
-323 0-0174 0-0396 0-1605 0-0174 0-0486 
-323 0-0285 0-0717 0:2460 0-0285 0-0950 
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SOME FOURIER TRANSFORMS IN PRIME-NUMBER 
THEORY 
By A. P. GUINAND (Ozford) 
[Received 15 August 1946] 

1. Introduction 

It has been shown that, if the Riemann hypothesis is true, then the 
relationship between the distribution of the non-trivial zeros of the 
Riemann zeta-function and that of the logarithms of the powers of 
prime numbers can be expressed in various ways closely connected 
with Fourier and Hankel transformations.+ 

In the present paper I[ discuss two pairs of Fourier cosine-transforms 
which illustrate one aspect of this relationship. I have previously 
discussed a related pair of Hankel transforms,{ but it is of some 
interest to construct a result which only involves the simpler Fourier 
cosine kernel. In each of the present pairs of transforms one function 
has simple discontinuities of magnitude (27)!/x when the argument x 
passes through a zero of {(}+-¢x), while the transform has simple 
discontinuities of magnitude 1/mp!” when the argument passes 
through a value of log p”, where p” is a positive integral power of a 
prime p. 

[ also derive an alternative proof of an infinite-series formula for 
N(T), the number of zeros of f(s) in 0 < I(s) < 7. 

All the results described above require the assumption of the 
Riemann hypothesis. Some simpler analogous results requiring no 
unproved hypothesis are given in the last section. 

2. First pair of transforms 


Suppose that the Riemann hypothesis is true, and let }+7y, 


(n 1, 2, 3,...; O < y, < Y,4,) run through the non-trivial zeros of 


vn 
f(s).§ Then 


x 


* AN N(x ee 
> i -{Se=--4 [vos (2.1) 
vaca?” : , a J 


1 
+ A. Wintner, Duke J. of Math. 10 (1943), 99-105 (99), and A. P. Guinand, 
Proc. Lond. Math. Soc. (to appear shortly), referred to in the sequel as (A). 
t (A), Theorem 1. 


§ If }+-7y, is a multiple zero of {(s) of order r+ 1 then we put 


Ya = FYaua = «> = Yaur 
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since 

Nowy N(x) = = ing—-——4-=-4- R(x), 
27 au 627. (8 

where R(x) = O(log x). 


Hence (2.1) becomes 


x 


t dt dt dt 
Hels { (wes, mit |? Wotets Nf 


1 1 seer 
= —log*x— 2 log 27 log x— _ (1++log 27)+ 


a7 


v 
= 
mil 


1 
5~ log 


where 


1 


Dae 


log 27 log a. 


: 1 1 
= lim ae = log*a+-5 
x00 Yn 7 
Yn 


_ 


Now put 


\ 


log 27 log ee 


] 


27 


ae ] = 
— —— log*x+ 
yn<u Yn ‘ 
where the dash indicates that the terms y,, = x, if they occur, are to 
be halved. Then 

Yn 


g F'(t)cos at dt 


+1 
Yn 


] 1 
1— +... — cos xt dt — 
Y2 Yn < 
Yn 


Yn 'N 
| log*t cos at dt +- log 27 2m | log ¢ cos at dt —2k | cos xt dt 
0 0 


0 


J 
— = 
2m J 


t E. C. Titchmarsh, The Zeta-function of Riemann (Cambridge, 1930), 4. 
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n 


1 7 ‘ 
$f... +2) (sin vy, ,,—sin xy,,)— 


Yy 
] 1 dt 
— —  log*y, sinryy+ — | logtsinat— 
5 ne Oey v— | g 7+ 
0 


bios ; bo fs ee ee 
+ — log 27 log yy sinxyy —— log 27 | sinat—— — sin 2yy 
7X 7X 3 t x 


0 


aS a | 1 
+3 —— — log*yy +— log 27 log yy -1| - 
Yn 4r 2a 


+ | logesinat— Tog 2x | sinat t+ o( 28%), (2.3) 
TX t 72x t YN 
0 0 


Now, by (2.2), the expression in the brackets { } is o(8%), Also 
YN 


ie] in8) 


ne) 
- 6 _ ae . ae 
logtsinat— = | logusinu——logz | sinu— 
t U u 
0 0 


0 
= —}n(C+log2), 
where C is Euler’s constant. Hence (2.3) becomes 
9 N 47 “ . 
—- 3 nocd ie (C+log 27a)+ o(?82) 
a o=1 Yn ot YN 
Now it is known thatt 
£’(s) Ss’ logp _y'- S co er 
ro § > = ee —” lS = D) = aa 
C(s) <P i—s 2, te ap? 
where p runs through the non-trivial zeros of {(s); s A 1, —2q, p; and 
~y > 1. If we put s = }, p= $+iy,, y = e*, then it follows, after 
some manipulation,{ that the series 
i 2) . 
sin xy, 
Yn 


2 
n=1 
+ E. C. Titchmarsh, The Zeta-function of Riemann (Cambridge, 1930), 81. 


t See (A) for details. 
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converges to the sum 


, log2 

ee a 4sinh $x+ } logeoth }a+ 

mlogp<zx p ; 
+arctan e-*—}3C—1n—}log8r. (2.5) 


Substituting (2.5) in (2.4) and making V tend to infinity, we have, 


after rear rangement, 


>a 


1 
(- . fr t)cos xt dt 


{ + logp 
| 


— 4sinh }7+-4lo 
yim ” ™ 


La 


anh jz 
g | + ({7—arctan em) 
4° 


ace 

= G(x), say. 
That is, G(x) is the Fourier cosine-transform of F(x), and the integral 
converges in the ordinary sense. 

Further, it follows from (2.2) that F(a) belongs to L7(0,00). Hence 
G(x) also belongs to L7(0,00), and, since F(x) is of bounded variation 
in any finite interval excluding the origin, it follows by a theorem 
of Titchmarshy that 


0 


r\1 


F(x) = -\ | G(t)cos xt dt. 


7 


0 
Thus we have 


THEOREM 1. Jf the Riemann hypothesis is true, and 


7 of 1 | 
F(x - (2 if i}, o2 =! o2rloga—k 
(x) = ( ” > : — log «+ — log 27 log 2 rf 


yn<a 


where k is chosen so that lim F(x) = 0, and 


ro 


+ logp tanh tx 


fe ae sts p—hx 
— + 47— arctan e . 


— 4sinh a+} log 
42 


4m 
the 
miog pcx I 


—>0O 


then, for x > 0, ; . | G(t)cos at dt, 
0 


oO 


P | F(t)cos at dt. 
0 


+ E. C. Titchmarsh, Fourier Integrals (Oxford, 1937), Theorem 58. 
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3. Second pair of transforms 


We use the following lemma:t 


Lemma. If f(x) and g(x) belong to L*(0,00) and are a pair of Fourier 
cosine-transforms, then the functions 


fo) zx 
dt 1 
fe) [ 1%, ae) | oe at 
x 0 
are also a pair of Fourier cosine-transforms belonging to L*(0, 00). 


Now suppose that {a,,} is an increasing sequence of positive numbers 
tending to infinity, that {a,} is another sequence not necessarily 
increasing or positive, that A(x) is everywhere differentiable, and that 


fee) = >| >" 4,—4e@) (3.1) 


On<zr 
belongs to L?(0,00). If we choose M so that ay_,; << # < ay, and put 
N > MU, ay = T, then 


T om 





yy at it 
fO>= (0, +g} + ya) | ot 
N-1 On+1 - 
dt dt 
+ az (ay $0,444) | 7 | AWS 
n=M F4 : 
M-1 N-1 
— a" 2 an+ (a,+a,+...+a,) pes” a 
1 “M Xn On+1 
' n=1 n=M 
A(t)|* a 
+P -f 407 
a 1 e. “= 
" 25 a2, 
n=M “n=1 n= 
A(T) Ale) [gird 
( x — 
mae fA 
+A) 20) | aro 


x 


+ This follows immediately from E. C. Titchmarsh, Fourier Integrals 
(Oxford, 1937), Theorem 69. 
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> a,—Ata} | 5 | aot} 
T 


+>" be w4|—7 be {doa} (3.2) | 
Pa tn <r | 


Now the first expression in (3.2) is equal to f(x). Hence, making 7’ 
tend to infinity, we have 


a f 04 ee a Pi fa 4K, 





where 
i 
K = lim [Pe ~fa ot|—- {> «,—a(0)}| 
T+ Xn t rs f 
Using a similar notation for 6,,, Ree , and B(x), and putting 
oe, dt | 

2) = B'(j——-L 


where L is a constant, then, if g(a) belongs to L°(0,00), an argument 
similar to the above shows that 


ae} f a t) dt =~ 3 b,—B(x )} 


0 ar 


If we now put B,, = y,, 6, = (27)?, L = k, and 


then ¢(x) is equal to F(x) of Theorem 1, and 





.s J (2zr)*/ ,, x Ss 
g(e)—~ f gyae= “nay F og +s), 8.8) 
0 
where N(x) = >’ 1 = HN (x—0)+- N(x+0)}. 


YnQr 
Further, if we put 
= mlog p, a, = log p/p*", A(x) = 4cosh 4a—¢(2) 


Xn 


ly 
log (“Se +427—arctan e-**+- 4e-#, (3.4) 





and d(x) = 


noe 
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then f(x) is equal to G(x) of Theorem 1, and, by the lemma, the 
Fourier cosine-transform of (3.3) is accordingly 


a [10% - _-- > >” spn? | sine t—x+ [ ood, 
; mp m t t 
mlogp<z i f 
(3.5) 


where 
- 


K = lim - a | sinh ye. 
T->c0 mp ™ t 
i 


mlogp<T 


4m 


7 Ss “EP — 4cosh }7| + v4 © or] 


ions 0 
milogp<T I 


1 
- f vot 2 | sinh t+ 
i 


0 


T 

1 dt 1{ 
a, | pest 
| sinha 77 


+ lim ; 
mp?” 


a, . log p<T 
[f we write / for the last limit above, then 


x= [oot 42 fannie 


Substituting this in (3. 5) we have 


4 a) 
Fo ry 1 . dt dt 
angie oil | a oo ‘eC. 
a | 10! mpi | me's [¢ @ t 
% < 0 z 
(3.6) 


milogp< 


Now the function | 0 (3.7) 


does not reduce to any simple expression, and it is more convenient 
to eliminate it. The Fourier cosine-transform of (3.7) is 


[ cosa ae{ $’(u) 
0 i 


ne) 


t=a 
2(2)ifsinae f wont] 42(2}* { sosinar, 


t 
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Differentiating (3.4) we find that 
1 2 
$'(z) = j(cosech Le — =)+ i sech }a—2e-#*. 


On substituting this in (3.8) the integrated terms vanish, and (3.8) 
becomes 


1 


1 /2\3 ~ 2 2\3 f _ dt 
—|-—}* ! cosech 1t— 2| sin ae ao “\? | sech Jtsinat— — 
42\7 _ 4x . t 
0 0 
2/2\4 
feo 
v 


= (4 I,+41,—2]3;), say. (3.9) 


ait 
noe 


Now,} if R(z) > —1, 





. __p-a 
log T(1-4 2) = | ig? ae e \ dt 


zeé 
| e—1 jt’ 














Putting 2 = —4-+ 2 and taking the difference we get 
1 YH F 
5; 08 (e+ a = amI(4+ iz) 
~T(4—12) . 
r ut 
= | { xe — ; sin zt! . 
J | e—] Jt 
Lee] 
1 2) = __ sin at| dt 
cosech ae in nt — 
a iJ” + Jl t Jt 
tte f fsint sina) a P fsint etl (3.10) 
. \ ¢ at | t wed 
0 
Now the second expression in (3.10) is a F seed integral, and is equal 
to xlogx. Further, for R(s) > —1, 
[ (@*-®sint—t*—e~) dt = —P(s—1)cos }s—T(s) 
0 


e+1)/, _ leon gon} 
1—s | 8 r 





+ This follows from Binet’s first integral formula for logIT(z). Cf. C. 
Stewart, Advanced Calculus (London, 1940), 493. 
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Letting s > 0, we find that 
sint ex e-| dt = Il, 
t t 
0 

and hence (3.10) gives 

I, = —2{am ['(4+ix)—axlogx+2}, (3.11) 
where am I'(}+-72x) is defined by putting am ['(}) = 0 and continuing 
analytically along the line $+-72. 

Nowt | sech }¢cos xt dt = msech zz. 
0 

Integrating with respect to x we have 


I. 


> 


, lt ‘ 
= [ sech sisin at = “| sech ru du = arctan(sinhzx). (3.12) 


0 0 


Further L= | e-¥sin nt = arctan 22. (3.13) 
0 


Substituting (3.11), (3.12), (3.13) in (3.9) we find that (3.7) and 
J] 


(27 


r [| —{am I'($+i2)—ax logx+2}+ } arctan(sinh 7x)—4 arctan 22] 
i i 
(3.14) 
are a pair of Fourier cosine-transforms of L*(0,00). Adding (3.7) to 
(3.6) and (3.14) to (3.3) we obtain another pair of Fourier cosine- 
transforms of L?(0,00). These functions are also of bounded variation 
in any finite interval, and hence it follows as in Theorem 1 that the 
Fourier integrals concerned converge in the ordinary sense. The 
result is: 
THEOREM 2. If the Riemann hypothesis is true, andt 
4 1 
’ 1 , t 
H(«) = —__—2 | sinh }¢— —I, 
mpi ie 
0 


milogp<x 


+ E. C. Titchmarsh, Fourier Integrals (Oxford, 1937), 177 (7.1.6). 
If necessary we may substitute 


id 
2 [ sinh "e = li(e#*)—li(e-*) 


0 
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where 
T ; 1 
: 1 , dt i og p : 
[= im ~—2 | sinh Wo —7| ~ —20 

iaeie m a mp? 0 t T -* p i 
and (Sie 
h t 
K(x) = 7) (Ny@)—amP(h+ iv) += log 29+ | 
z | 2a 
1 , 2 | 
+ — arctan(sinh 7x) —— arctan 22], 

47 7 
9\1 > 
then, for x > 0, H(z) = (=\* | K(é)cos at dt, ; 
7 t 
0 ’ 
0 
9\1 } 
and K(x) = (=) | H(t)cos xt dt. (3.17) | 
7 i 
‘ 


4. The formula for N(z) 
As before, putting «, = mlogp, ay = T,a, = log p/p, the right- 
hand side of (3.17) is 


ae) t 

| - Gn _9 | sinh du -| cos at dt 
: 

0 0 


mt 


. N-1 On+1 
— (7\3 3; > ot “ [ sat dt — 
(=) sim (@ + te et. ; cos xt d. 


n=1 


Xn 


d t 
i . ae oe 
un G | cos zt dt | sinh 4a ——-sinaT 
3 og” BD 


0 0 


N- 
2\3 , [1 Ay | I Ly ; 
= {-}° lim 15 (2 es t.. Pe (sin xa,,,;—SiN Ya,) — 
TT} No] & Xn 
n=1 
. 
9 
2 ; ‘ ; du tl. 
—2 f sinh (sine? —sin aw) Mt —Lsine| 
x © 2 


0 


1/2\3 .. ae @ 
= —{—]° lim — > einza,+ 
X\T;) Now Xn 
Pom 4 


pe r r 
a, du du} 
—*—2 | sinh }u << +2 | sinh usin zu— 

at Oy u u 


+ sin ap >=- 


1 0 
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By (3.16) this is equal to 


l 2\b a, du 
=(=) ain »| - > i. sin va, + foram 


. N 
sn) > a,—20i"| 2 2\ ao 
T | x . 


n=1 





y 


T 
1/2\3 ,. ] 

— —-—(—}* lim > sin(am log p) — once 
xu\n) T-+0 mp u 

mlogp<T 0 
sinaT’ ] wd 1/2\3 
ais a | “Ef — 2Qe47 ] —+(} arctan 22, 
mlogp<T Pp om 
und by (3.17) this expression is equal to K(x). Substituting the value 
’ K(x) from Theorem 2 and rearranging the terms, we obtain the 





result: 
THEOREM 3.} If the Riemann hypothesis is true and x > 0, then 


1{N(a—0)+N(a-+0)}— (= log 5 -5 





r — du 
a im| > a sin(em log p)— | ei” sin 2u— — 
7i To Pe : Uu 
sinxT’ log p , 1 ‘ 
—_— — — 2a +5—{am I'(}+ix)—2 log x+2}— 
mlogp<T r eal 


1 F 1 
— — arctan(sinh 72) + — arctan 22, 
4ar 7 


where amI(4+izx) is defined by taking amI(}) = 0 and continuing 


['(s) along the line s = $+-12. 
Theorem 8 is, in a sense, analogous to (2.5), and the argument of 
this section can be reversed to deduce Theorem 2 from Theorem 3. 


5. Simpler pairs of transforms 


The pairs of transforms discussed in §§ 2, 3 have simpler analogues 
with regularly spaced discontinuities. For example, the functions 


A>, al, o-| > 5 hes] (5.1) 


+ See (A), Theorem 2, for an alternative proof and discussion of the result. 
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are a pair of transforms with respect to the kernel 2 cos 27x. Further, 
if y(n) is a real primitive character modulo « (« > 1) and 


x(n) _ 9g, 
n 
n=1 
then the functions 

: l~ x(n) 
- n), sS— > *— 
3 x(n) é n 

NSE oe 


are a pair of transforms with respect to the kernel 2«~! cos(272/«) if 
x(—1) = 1, or with respect to the kernel 2«~! sin(27x/«) if 
x(—1) = —1. 

These results are easily proved by the method of § 2, using an 
ordinary Fourier series in the place of (2.5). The pair of transforms 
(5.1) can also be derived as a limiting case of an earlier result. 

t A. P. Guinand, J. of London Math. Soc. 14 (1939), 97-100. Let s > 1—0 
in (1). 
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to apply statistical methods in various fields. . . . Such must 
welcome Mr. Kendall’s book warmly.’—Nature. 
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